arXiv:1503.00476v2 [cond-mat.soft] 25 Sep 2015 
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We demonstrate that the formation of bicontinuous emulsions stabilized by interfacial particles 
(bijels) is more robust when nanoparticles rather than microparticles are used. Emulsification via 
spinodal demixing in the presence of nearly neutrally wetting particles is induced by rapid heating. 

Using confocal microscopy, we show that nanospheres allow successful bijel formation at heating 
rates two orders of magnitude slower than is possible with microspheres. In order to explain our 
results, we introduce the concept of mechanical leeway i.e. nanoparticles benefit from a smaller 
driving force towards disruptive curvature. Finally, we suggest that leeway mechanisms may benefit 
any formulation in which challenges arise due to tight restrictions on a pivotal parameter, but where 
the restrictions can be relaxed by rationally changing the value of a more accessible parameter. 

PACS numbers: 81.16.Dn, 82.70.Dd, 64.75.Xc, 68.08.Bc 


I. INTRODUCTION 

The directed assembly of colloidal particles enables 
the design of novel soft materials with bespoke 3D ar¬ 
chitectures. The desired assembly route can be selected 
by adjusting the interparticle interactions. For example, 
the electrostatic interaction between oppositely charged 
particles can be tuned to obtain ionic colloidal crystals 
rather than irreversible aggregation [l|. An alternative 
approach employs templates to guide particle assembly 
towards a target structure. For instance, sedimenta¬ 
tion of microparticles onto structured solid templates has 
been used to direct colloidal-crystal assembly Q and bi¬ 
nary crystals of nanoparticles have been grown via liquid- 
air inter facial assembly Q- In both cases, the interac¬ 
tion between the assembling particles and the template 
is crucial: pattern-lattice mismatches of ~ 10% already 
cause crystal defects and liquid subphase properties sig¬ 
nificantly affect crystal quality [2], Q. 

A startling case of liquid templating is the formation 
of bicontinuous Pickering emulsions Q, i.e. bicontinuous 
interfacially jammed emulsion gels or bijels (Fig. Oja)) 
0 H , which have been suggested for applications in fuel 
cells, microfluidics and tissue engineering [Ml. Bijel 
formation typically proceeds via spinodal demixing of a 
binary liquid containing colloidal particles (Fig. [U(b)), 
which can arrest the phase separation by forming a 
jammed monolayer at the liquid-liquid interface. As in 
the cases discussed above, template-particle interactions 
are essential: bijels are only formed if the particles are 
(almost) neutrally wetting, otherwise emulsion droplets 
are formed [9]. The parameter that quantifies this inter¬ 
action is the contact angle 6 , which is a measure of the 
particle’s position relative to the liquid interface: 0 = 90° 
is neutral wetting (Fig. [U(c)). Unfortunately, tuning the 
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mean value of 6 is non-trivial and restraining its variance 
is harder still, making bijel formation challenging. 

Ostensibly, reducing particle size r given a fixed final 
bijel-channel width Lf r (Fig. DJa)) would only make 
matters worse, as scaling down r in a close-packed mono- 
layer of particles with fixed 6 % 90° requires a commen¬ 
surate reduction in the local radius of curvature of the 
interface [15]. In other words, for a given non-neutrality, 
one might expect smaller particles to locally demand a 
more strongly curved interface and hence be more disrup¬ 
tive to bicontinuity on a chosen scale Lf >> r. However, 
this ignores the particle-size dependence of the stiffness 
of the particle-laden liquid interface, which might specif¬ 
ically aid small particles in overcoming off-neutral wet¬ 
ting. 

In this paper, we experimentally explore the effect 
of particle size on bijel formation. We find that bi¬ 
jels are formed more robustly when nanoparticles rather 
than microparticles are used: nanospheres allow mini¬ 
mum heating rates two orders of magnitude slower than 
microspheres, with the latter stabilizing droplet emul¬ 
sions rather than bijels at slow rates. We discuss our 
results in the context of mechanical leeway, i.e. inter fa¬ 
cial particles that are smaller lead to a less rigid inter¬ 
face between the two liquid phases, resulting in a smaller 
driving force towards disruptive curvature. Finally, we 
discuss the implications of leeway mechanisms in the (di¬ 
rected) self-assembly of functional formulations based on 
particle-template or even particle-particle interactions. 


II. MATERIALS AND METHODS 

A. Materials 

For particle synthesis, tetraethyl orthosilicate (TEOS, 
> 99%, Aldrich), 35% ammonia solution (reagent 

grade, Fisher Scientific), ethanol absolute (VWR Chem¬ 
icals), fluorescein isothiocynate (FITC, 90% isomer 1, 
Aldrich) and (3-aminopropyl)triethoxysilane (APTES, 
99%, Aldrich) were used as received. For bijel prepa- 
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FIG. 1. (color online), (a) Schematic of a 2D slice through 
a 3D bijel: two continuous liquid channels (A black and B 
magenta), stabilized by a jammed layer of interfacial particles 
(yellow). Lf: final channel width, (b) Coexistence curve for 
the water-lutidine (W-L) system (CP: critical point) [T3 ]. 
Vertical arrow: bijel formation, i.e. a homogeneous mixture 
of W-L at the critical weight fraction (X) is heated from room 
temperature to 45 °C (or 50 °C). Spinodal demixing results in 
two phases (A/B), with compositions given by the horizontal 
tie-lines, (c) Schematic of the contact angle 6 (for 0 — 90°). 


ration, 2,6-lutidine (> 99%, Aldrich) and Nile Red 
(Aldrich) were used as received; distilled water was run 
through a Milli-Q (Millipore) filtration system to perform 
deionization (to a resistivity of at least 12 MD • cm). 

Here, we formed (bicontinuous) Pickering emulsions 
by spinodal demixing of the binary liquid water-lutidine, 
heated at various rates in the presence of colloidal parti¬ 
cles. Note that the water-lutidine (W-L) interfacial ten¬ 
sion 7 wl is temperature-dependent and orders of mag¬ 
nitude lower than that of typical water-alkane systems. 
According to Ref. Q , 7 wl ranges from ^0.01 mN m 1 
at 34.2 °C (just above the lower critical solution temper¬ 
ature of 34.1 °C) to ~ 0.4 mN m -1 at 46.0 °C. During 
slow heating at 1 °Cmin _1 , it takes about 6 s to get from 
34.1 °C to 34.2 °C and about 12 min to get to 46.0 °C. 



(b) w 



• • 


FIG. 2. Transmission Electron Microscopy (TEM) images 
of the Stober silica (a) microparticles and (b) nanoparticles 
used in this study. The particles have an average radius of (a) 
348 nm and (b) 63 nm; polydispersity (a) 6.4% and (b) 15% 
(analyzed images taken at different magnification to improve 
statistics). Scale bar: (a) 1 /xm and (b) 200 nm. 


B. Particle Synthesis 


The particles used in this study were synthesized using 
the Stober method 16], modified to include the dye FITC 
via the linking molecule APTES [17]. For the micropar¬ 
ticles (MPs), a dye mixture of 0.584 g APTES, 0.107 
g FITC and 4.0 mL ethanol was prepared overnight by 
stirring. The following day, a reaction mixture of 1.5 
L ethanol, 186 mL 35% ammonia solution and 60 mL 
TEOS was prepared, and the dye mixture added. The 
entire reaction mixture was kept in a refrigerator for 24 
hours at ^ 10 °C. This resulted in particles with a ra¬ 
dius of 0.36 fi m as measured by Dynamic Light Scattering 
(DLS) and 0.35 /im according to Transmission Electron 
Microscopy (TEM, Fig. [2ja)). 

The nanoparticles (NPs) were synthesized in a sim¬ 
ilar fashion to the MPs, except the reaction tempera¬ 
ture was 25 °C and the concentration of dye mixture 
was increased to take account of the increase in surface- 
to-volume ratio which accompanies a decrease in particle 
radius. This is an important consideration, as it has been 
shown that the presence of APTES on the silica surface 
is crucial for meeting the neutral-wetting requirement in 
the W-L system mui - it has been suggested that 
the surface decorations act to disrupt the wetting layer 
of lutidine which spontaneously forms around the parti¬ 
cles when approaching the phase separation temperature 
[2()1 . 2l]]. For the NPs, DLS returned a particle radius of 
0.08 fi m and TEM returned 0.063 /im with a polydisper¬ 
sity of 15% (Fig.[2fb)). We have confirmed that the NPs 
(1.51 ±0.06 g mL -1 ) have a lower density than the MP 
(1.63 ± 0.03 g mL -1 ) (density meter, Anton Paar, DMA 
4500), presumably due to the higher dye concentration 
by volume j^|, which could lead to enhanced shrinkage 
in the vacuum of the TEM [23]. The 4.5 x decrease in 
DLS particle size closely matches the 4.5 x increase in 
APTES concentration compared to the MP synthesis, so 
the NPs and MPs are expected to have identical surface 
chemistries. 
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To remove excess APTES and FITC from the synthe¬ 
sis product, the particles were washed by repeated cen¬ 
trifugation/redispersion: 2x ethanol then 4x water for 
the MPs and 2x ethanol then 8x water for the NPs. 
Subsequently, the particles were pre-dried at room tem¬ 
perature in a fume hood and ground with a mortar and 
pestle. Prior to sample preparation, particles were dried 
at 20 mBar and 170 °C (no more than 100 mg per vial 
and no more than 3 vials at the same time) [jjj. This 
removes surface-bound water and may cause moderate 
dehydroxylation of the silica surface (24|. The drying 
time was tuned to optimize bijel quality as assessed by 
visual inspection of confocal micrographs; dried particles 
were stored in a desiccator in the presence of a silica gel. 


C. Sample Preparation 

First, dried particles were dispersed in deionized wa¬ 
ter by ultrasonication (Sonics VibraCell). The MPs were 
sonicated for (2 x 2) minutes at 8 W with (2 x 10) s of 
vortex mixing in between. To ensure proper redisper¬ 
sion, NPs were additionally sonicated for (lx 10) min¬ 
utes at 8 W and vortex mixed for (1 x 10) s. Lutidine 
was then added to give a mixture with a critical composi¬ 
tion, i.e. a mass ratio of W:L = 72:28 (Fig. Hfb)) [llj, so 
that spinodal decomposition would be (at least initially) 
the preferred phase separation mechanism. To allow con- 
focal imaging of the lutidine-rich phase, the fluorescent 
dye Nile Red had been added to the lutidine at a con¬ 
centration of around 10 fi M (we checked that Nile Red 
partitions into the lutidine-rich phase and that concen¬ 
trations as low as 1 g M gave similar bijels). The sample 
mixture was transferred to a glass cuvette (Starna 21- 
G-l with pathlength 1 mm) and placed inside a metal 
block, which was itself placed inside a temperature stage 
(Instec, TSA02i). Emulsification via liquid-liquid demix¬ 
ing was initiated by heating the sample to a target tem¬ 
perature above the lower critical solution temperature 
(LCST) of 34 °C. 

Slow heating (< 5 °C min -1 ) was achieved by pro¬ 
gramming the temperature stage to ramp the temper¬ 
ature T at the desired rate T, from room temperature 
(« 20 °C) to 45 °C. Heating rates were extracted 
from the (T, t)-graphs produced by the stage software 
and we have used a thermocouple to ascertain that at 
these slow rates the sample temperature does not lag 
the stage temperature; estimate of corresponding error 
in heating rate <jf = 0.1 °C min -1 . For a heating rate 
of 17 °C min -1 , we adopted a method from Ref. the 
temperature stage and metal block were pre-warmed to 
45 °C or 50 °C and the room-temperature cuvette was 
inserted. We have confirmed this heating rate by mea¬ 
suring the time it took to reach phase separation at the 
LCST of 34 °C from room temperature; estimate of cor¬ 
responding error in heating rate Of = 3 °C min -1 . For 
higher heating rates, the cuvette was placed on top of 
a small cardboard box (to prevent thermal conduction 


away from the cuvette) inside a microwave (DeLonghi, 
P80D20EL-T5A/H, 800 W, set to “auto-defrost 100 g” 
i.e. 40%) 0. The sample was irradiated for 5 s (or 
6 s) and then quickly transferred to the temperature 
stage at 50 °C. We have checked by visual inspection 
that the sample remained opaque (i.e. phase separated) 
upon transfer from the microwave to the temperature 
stage. The corresponding heating rate was calculated as 
(50 °C — 20 °C) /5 s = 360 °C min -1 , with an estimated 
error of 30 °C min -1 . 


D. Characterization &; Image Analysis 

During or after emulsification, samples were imaged 
using fluorescence confocal microscopy. Fluorescence ex¬ 
citation was provided by a 488 nm laser (for FITC) and 
a 555 nm laser (for Nile Red); emission filters were used 
as appropriate. The two liquid domains could be dis¬ 
tinguished by detecting the fluorescence of the Nile Red, 
while the location of the particles could be determined 
by detecting the fluorescence of the FITC. 

To extract the bijel channel width L from 2D confocal 
microscopy images, a pixel-based correlation function al¬ 
gorithm was run on the Nile-Red channel using Matlab. 
The algorithm constructs a radial distribution function 
g (r) by multiplying pairs of pixel intensities, plotting the 
values against the distance between the pixels, and then 
taking an average; the bijel channel width or character¬ 
istic length scale is then taken to be the location of the 
first minimum in the plotted g (r ) jj25j. For the final bijel- 
channel width Lf, this process was repeated on at least 
three separate images of the same bijel sample and an 
average was taken. The standard deviation of measure¬ 
ments made on several images of the same sample was 
taken as the error gl ~ 3 gm. 


III. RESULTS 

We begin by comparing (bicontinuous) Pickering emul¬ 
sions formed by spinodal decomposition of W-L mixtures, 
containing either nanoparticles (NPs) or microparticles 
(MPs), upon heating at various rates (Fig. [TJb)). Fig. [3] 
presents a confocal-microscopy overview of the structures 
obtained for two different particle radii and three dif¬ 
ferent heating rates. In all panels, the fluorescently la¬ 
beled particles (yellow) appear at the liquid-liquid inter¬ 
face between the water-rich phase (black) and the fluo- 
rescently labeled lutidine-rich phase (magenta). Samples 
prepared with MPs show bicontinuous structures only for 
fast heating (Fig. [3fa)), whereas slow heating results in 
discrete droplets (Fig. [3{c)). In contrast, NPs invariably 
yield a percolating interface with both signs of curva¬ 
ture (Fig. (3jd-f)), which is an imperative characteristic 
of a bijel; note that slow heating with NPs (Fig. [3^e,f)) 
seems to yield a relatively higher number of thin necks 
compared to fast heating (Fig. [3fd)). 
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FIG. 3. (color online). Fluorescence confocal micrographs of 
final-state emulsions of water and lutidine (magenta), formed 
using various heating rates (dT/dt), stabilized by (nearly) 
neutrally wetting particles (yellow) of radius r. Particle vol¬ 
ume fraction is (a) 2.6%, (b-c) 2.2% and (d-f) 0.7%. Es¬ 
timated relative error in heating rate erf < 18% (Sec. HTCl) . 
Scale bars: 100 /im. See Appendix lAl for sample homogeneity. 



Next, we compare the kinetics of bijel formation using 
MPs vs NPs, to explain the discrepancy in the struc¬ 
tures obtained after slow heating (Fig. [3]). Fig. 0] shows 
selected confocal micrographs from time-series recorded 
during slow heating in the presence of MPs vs NPs. Using 
MPs (Fig.0]^a-d)), the interconnected domains present at 
t = 2 s have pinched off by t = 4 s, resulting eventually in 
particle-stabilized droplets. By contrast, when using NPs 
(Fig. 0te-h)), connectivity is maintained until the struc¬ 
ture is arrested, resulting in a bijel. Though we observe 
thinning of necks, we cannot find a convincing pinch-off 
event between Fig.0|(g) and Fig.0^h). Note that we have 
also observed droplet formation via secondary phase sep¬ 
aration (Appendix [B]), but this does not seem to be a 
pivotal effect, i.e. it can both happen and fail to hap¬ 
pen irrespective of bijel formation failing or succeeding 
[6. [l9|, 0, [27]. This suggests that MPs fail to produce 
bijels via slow heating, because depercolation via pinch- 


FIG. 4. Time sequences of confocal fluorescence micrographs 
showing water-lutidine mixtures containing (nearly) neutrally 
wetting particles of radius r (white) during slow heating 
(1 °Cmin -1 ). Particle volume fraction is (a-d) 2.1% and (e- 
h) 1.8%. Note (c,d) the depercolation via (encircled) pinch-off 
events and (e-h) the formation of a bijel (also verified down 
to (^np = 0.7% (Appendix [C])). Scale bars: 100 /xm. 


off events occurs before the interfacial particles jam and 
lock-in the bicontinuous structure. 

To quantify the coarsening observed in Fig. 0J we used 
image analysis to extract the channel width L (Sec. mm. 
Fig. [5fa) shows that the coarsening in the presence of 
MPs is similar to coarsening without particles, until 
t = 5 s when the bicontinuous structure has failed and 
MP-stabilized droplets have appeared. Coarsening in the 
presence of NPs initially follows the behavior of the W-L 
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FIG. 5. (a) Measured channel width L vs time t during spin- 
odal demixing upon heating at 1 °Cmin -1 of a critical mixture 
of water-lutidine without (W-L) and with (MPs) 0.36 /im ra¬ 
dius microparticles or (NPs) 0.08 /xm radius nanoparticles, the 
latter resulting in a bijel. (b) Corresponding coarsening speed 
vl for the NP data. The dashed vertical lines enclose the jam¬ 
ming time Atj «4s. Estimated error in L is cfl ~ 3 /im. 


mixture without particles, but then levels off. As bijel 
formation at 1 °C min -1 fails with MPs (and without 
particles), Fig. Efb) only shows the coarsening speed 


vl 


Lti ~ itj-i 
ti ti —1 


(i) 


in the case of NPs; note that vl goes through a maxi¬ 
mum at t « 3 s and is (more or less) 0 after t ~ 7 s. As 
discussed below, we refer to the time between the max¬ 
imum in vl and its levelling off as the ‘jamming time’ 
-Vi- 


IV. DISCUSSION 

Having presented our experimental results, we first dis¬ 
cuss how bijel formation can fail and how particles with 
off-neutral wetting can promote bijel failure. Simula¬ 
tions of spinodal demixing without particles in 3D, in the 
viscous hydrodynamic (VH) regime relevant here, have 
shown that depercolation proceeds via thinning of liq¬ 
uid channels followed by pinch-off events [28|. Neutrally 


wetting particles can halt the demixing by attaching to 
and jamming at the liquid interface [HfE, 9]. However, 
off-neutral particles induce a spontaneous curvature Co 
when attached to liquid interfaces [15, 19, 29, 30]. This is 
because they are pushed together as coarsening decreases 
the interfacial area, while the interparticle contacts are 
not situated at the liquid interface (where they would be 
for 6 = 90°). As bijels have empirically been shown to 
feature average mean curvature (H) = 0 [10, EH, an y 
Co 7^ 0 is expected to disrupt bijel formation. 

Note that secondary nucleation, i.e. the formation of 
new droplets during spinodal decomposition, was not ob¬ 
served in the above-mentioned simulations jgl, [H, 3(j, 
presumably because the quench was instantaneous 27| . 
Secondary nucleation during bijel formation has previ¬ 
ously been observed in experiments and attributed to the 
finite rate of temperature change BM- However, it has 
not been suggested that secondary nucleation is respon¬ 
sible for bijel failure, rather it results in droplets inside 
bijel channels |6] or even droplet-reinforced channels [27j | . 

Intriguingly, our results show that bijel formation fails 
during slow heating with MPs, whereas it succeeds with 
NPs that were designed to have similar wetting. The NP 
contact angle could simply be closer to 90°. However, this 
does not agree with our observation that NPs allow bijel 
formation over a wider range of drying times, which is 
expected to correspond to a wider range of contact angles 
jl9| . Our fluorescence confocal time series suggest that 
MP bijels fail due to depercolation via pinch-off events. 
Pinch-off events may also occur for NPs: they can even be 
observed in 3D simulations of successful bijel formation 
[Hf. However, we suggest that NPs sufficiently suppress 
the number of pinch-off events to allow successful bijel 
formation. 

In order to explain why NPs facilitate bijel formation, 
we have found it particularly illuminating to consider the 
particle-size dependence of the “driving force” F towards 
Co (Appendix [Dj i.e. away from (H) = 0 for bijels 
l3I| . The bending-energy density of the particle-laden 
interface is 


w = 2k {H - Co) 2 , (2) 

where k is the effective bending modulus of the interface 

mm, so 


F = — = -4k(C 0 -H). (3) 

Dimensional analysis suggests that Co oc — 1/r and 
k oc 7wiT 2 , which is backed by analytical calculations 
for spheres on a spherical cap [15|. As here Lf >> r, and 
so \H\ ~ 1/Lf < |Co|, we approximate Eq. (j3]) as 

F « -4 kCo oc ywiT • (4) 

Thus, NPs demand a more strongly curved interface 
(Co oc —1/r), but the driving force towards that cur¬ 
vature is smaller (F oc r). 
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To assess to what degree a smaller driving force can 
facilitate bijel formation, we compare the disruption time 
Atd to the jamming time Atj; bijel formation can succeed 
if Atj < Atd- For the NPs, we can estimate the jamming 
time from Fig. [5jb) . We define the jamming time as 
Atj = tf — ti n , where ti n is the time at which the jamming 
starts causing a decrease in the coarsening speed vl - the 
peak in Fig. [5fb) - and tf is the time just before vl drops 
to zero. This gives Atj ? NP ~ 4s at a heating rate of 
1 °C min -1 . 

We cannot obtain the MP jamming time directly, since 
MP bijels fail at 1 °C min -1 . However, we expect the 
jamming dynamics to be dominated by the instantaneous 
area fraction of inter facial particles, which is independent 
of particle radius, as long as the final lengthscale is fixed 
( Appendix ID]) . As Lf (MP) > Lf (NP) in Fig. 01 i.e. <^np 
is 1.8% vs 0.7% in Fig. [3] and Lf oc r/(j) [l2j, we expect 
Atj ? MP > Atj ? np ~ 4 s. 

Conversely, we can estimate the MP disruption time, 
Atd,MP ~ 2 s from the time of occurrence of pinch-off 
in confocal images (Fig. 0(b-d)), whereas we cannot es¬ 
timate Atd,NP because bijel formation succeeds here for 
NPs. However, we can predict the scaling of At a with 
particle radius by balancing the driving (Eq. (j4|) and 
viscous-drag forces, to give 

Aid oc , (5) 

7wl r 

where r] is a bulk fluid viscosity and A >> r is the typical 
length scale of the disruption, which is independent of 
particle radius (Appendix ID]) . Given that the MPs are 
4.5 x larger than the NPs, and assuming effects of particle 
polydispersity and roughness are negligible, the inverse 
scaling of Ata with radius implies Atd,NP ~ 9 s. 

These time-scale estimates help to explain the ob¬ 
served patterns of bijel failure, i.e. they explain why 
np < Atd,NP but Atj ? MP > Atd,MP- To account for 
any possible dependence of Ata on the (final) channel 
width, we have also verified that bijel formation is suc¬ 
cessful with NPs at 1 °C min -1 for similar Lf, i.e. for 
0 np = 0.7 vol-% (Appendix 1C]) . It is worth noting here 
that, based on the scaling proposed in Eq. (J5j), we had 
expected that bijel formation would succeed with MPs 
at 5 °C min -1 . This is because it succeeds for similar 
Li with NPs at 1 °C min -1 (and the NPs are about 5x 
smaller than the MPs). As bijel formation with MPs is 
only barely successful at a higher rate of 17 °C min -1 , 
this suggests that an additional mechanism might be at 
play here; currently planned simulations and experiments 
may be able to address this in the future. Having said 
that, the mechanical-leeway mechanism proposed here 
does point in the right direction, i.e. it can explain why 
bijel formation is more robust when using NPs rather 
than MPs. 

As shown above, slow heating increases the importance 
of bypassing droplet formation. We have suggested here 
that NPs succeed in this because of their larger mechan¬ 
ical leeway, whereas MPs may fail under similar condi¬ 


tions (Fig. [3]). This also has technological relevance, since 
fast and homogeneous heating is challenging to achieve, 
putting severe restrictions on the choice of sample geom¬ 
etry and starting materials 9 }. Therefore, reducing par¬ 
ticle size could greatly facilitate formulation, especially 
when tuning particle surface chemistry is non-trivial (as 
is often the case), even though a naive expectation based 
on an optimal (static) wetting geometry would suggest 
exactly the opposite trend. 

This mechanical-leeway mechanism not only applies to 
bijels but to any liquid template for solid particles. More 
broadly, leeway mechanisms may well aid any formula¬ 
tion where challenges arise due to tight restrictions on 
a pivotal parameter, but where the restrictions can be 
relaxed by changing a more accessible parameter (here: 
particle size). This has important implications for the de¬ 
velopment of fabrication routes for advanced functional 
materials based on external templates. Moreover, it is 
potentially relevant to the design of any soft material 
with a bespoke architecture by adjusting particle interac¬ 
tions, e.g. crystallization of spheres with a size variation 
above the hard-sphere crystallization threshold (^ 12%) 
is possible by changing the ionic strength of the 
suspending medium so that the interparticle-interaction 
range is large enough [35). 

V. CONCLUSIONS 

We have shown that the formation of bicontinuous 
Pickering emulsions (bijels) via liquid-liquid demixing is 
more robust with nanoparticles than with microparticles: 
a wider range of heating rates can be used. In addition, 
our results suggest that bijel formation using microparti¬ 
cles fails at low rates because the bicontinuous structure 
decays into discrete droplets via pinch-off events. To ex¬ 
plain our observations, we have argued that inter facial 
microparticles with off-neutral wetting induce disruptive 
curvature, while nanoparticles of similar wetting benefit 
from a mechanical-leeway mechanism. In short, smaller 
particles give a smaller driving force towards disruptive 
curvature. 


ACKNOWLEDGMENTS 

M.R. is grateful to EPSRC for funding his PhD 
studentship. M.E.C. is funded by the Royal Soci¬ 
ety. J.H.J.T. acknowledges the Royal Society of Edin- 
burgh/BP Trust Personal Research Fellowship for fund¬ 
ing and The University of Edinburgh for awarding a 
Chancellor’s Fellowship. The authors are also grateful 
for financial support from EPSRC EP/J007404. Thanks 
to Paul Clegg, Michiel Hermes and Alexander Morozov 
for useful discussions. 



7 


Appendix A: Sample homogeneity 


Appendix B: Secondary nucleation 


In this Appendix, we present several fluorescence con- Below are confocal micrographs, corresponding to 
focal micrographs of a MP and a NP stabilized bijel, to Fig. 0J to illustrate secondary nucleation. 
demonstrate sample homogeneity. 



r = 0.36 urn r = 0.08 




FIG. 7. Confocal fluorescence micrographs from two time- 
series showing water-lutidine(white) mixtures containing 
(nearly) neutrally wetting particles of radius r during slow 
heating (1 °Cmin -1 ). Particle volume fraction is (a,b) 2.1% 
and (c,d) 1.8%. Note that droplets have appeared, presum¬ 
ably due to secondary nucleation, which has previously been 
observed during slow quenches [6, [27|. Scale bars: 100 /im. 


FIG. 6. (color online). Fluorescence confocal microscopy on 
bijels of water and lutidine (magenta), stabilized by nearly 
neutrally wetting particles (yellow), formed using microwave 
heating. Micrographs of a (a-c) microparticle (MP) and (d-f) 
nanoparticle (NP) stabilized bijel at three different positions 
(randomly chosen). Particle volume fraction is (a-c) 2.6% 
and (d-f) 0.7%. Scale bars: 100 /im. 
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Appendix C: Slow heating using NPs at 0.7 vol-% 


Appendix D: Timescales 


Here, we present a confocal micrograph demonstrating 
successful bijel formation at a heating rate of 1 °Cmin _1 
and a nanoparticle volume fraction of 0.7% (see caption 
of Fig. |4] and discussion after Eq. (j5j)). 


In this Appendix, we obtain approximate scaling re¬ 
lationships for the timescales of jamming and disruption 
during bijel formation. 



FIG. 8. Fluorescence confocal micrograph of a water-lutidine 
bijel, formed using a heating rate of 1 °Cmin _1 , stabilized by 
nanoparticles of radius 0.08 /im (white) at a volume fraction 
of 0.7%. Scale bar: 100/xm. 


1. Disruption Time 

Following Canham and Helfrich [32, [33], we start with 
the bending-energy density of a membrane 

w = 2 k (H - Co) 2 + k g K , (Dl) 

in which n is the bending modulus, H the mean cur¬ 
vature, Co the spontaneous curvature, kq the Gaussian 
bending modulus and K the Gaussian curvature. As¬ 
suming the topology of the surface does not change sub¬ 
stantially during the crucial stages of bijel formation, we 
omit the K term [H|: 

w = 2k{H-C 0 ) 2 . (D2) 

Next, we consider the (generalized) driving force F to¬ 
wards spontaneous curvature. Taking H as constant over 
a small membrane patch, 


77> _ 8W 

* ~ dH 


d 

dH 


2k (H - CoY 
= -4k (Co -H) . 


(D3) 


Eq. m resembles Hooke’s law for a spring with spring 
constant k = 4ft and extension u = (Co — H). The equi¬ 
librium position of the spring is H = Co, which is a min¬ 
imum as (d 2 w /<9i7 2 ) = An (which is positive for n > 0). 
Note that it has been shown empiricallythat the average 
mean curvature (H) = 0 for bijels jlQl. I31j . 

In order to understand how the driving force F scales 
with particle size r, we first consider how the spontaneous 
curvature Co and the bending modulus n scale with r. 
Co has units of inverse length (m -1 ) and is expected to 
scale as — 1/r, which is backed up by analytical calcula¬ 
tions for spherical particles on a spherical cap [l5|. In 
that geometry, the result can also be explained using a 
scaling argument: to keep the angles fixed, including the 
particle’s contact angle, both r and the radius of curva¬ 
ture R c of the spherical cap have to be reduced by the 
same factor, showing that 

c “ = • < D4 > 

Note that Co also depends on the particle’s contact angle 
0 and that Co = 0 for neutrally wetting particles (0 = 

90°) m. 

The bending modulus n has units of energy (J). As it 
is expected to depend on the W-L interfacial tension 7 wl 
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(units N m 1 ) and on the presence of the particles, one 
might guess 

ft oc 7 wl ^ 2 • (D5) 


This claim is backed up by analytical calculations of ft 
for a close-packed monolayer of spherical particles on a 
spherical cap [15]. 

In our experiments, the final bijel-channel width Lf >> 
r, so \H\ ~ 1/Lf |Co| (Eq. (jD4|) ). Combined with 
Eqs. droll and dEU), this means the driving force F scales 
with r: 

F = —4ft (C 0 — H) 


~ —4ftCo 
oc - 7 WL r 2 ■ 


(D 6 ) 


oc 7 WL r . 


In words, for the same binary liquid (twl) and a given 
off-neutral wetting (6 ^ 90°), the driving force towards 
the spontaneous curvature is smaller for NPs than it is 
for MPs, which can help explain why fabricating bijels 
is possible over a larger range of heating rates with NPs 
than with MPs. 

To gain a simple estimate of the disruption time A t d , 
which is the time it takes for the driving force F to cause 
so much curvature that bijel formation fails, we balance 
F with a viscous drag force: 

F = Tdrag 

oc rj\v , ( D7 ) 


As the Reynolds number Re <C 1 here, even when con¬ 
sidering motion at the scale of the channel width L, we 
can ignore the inertial term [37J: 

flit + 4/vn = 0 . (Dll) 

Re-writing Eq. (1D11[) results in an expression for the rate 
of change of curvature ( dH/dt ) 

u = 

d(Cp-H) _ 
dt 

d H _ 
dt 

Let us denote the time when the interfacial particles 
start interacting as t[ n . As at that time the bijel channel 
width L r, we can write 


_ 4 KU 

ft 

_ 4*(Co-g) ( D12 ) 

4 K (C 0 -g) 


(dH\ kC 0 


(D13) 


For bijel disruption to occur, the curvature H has to 
change by a threshold amount AH d ~ A -1 . For the 
disruption time, we can then write 


A -a Aft/d 

^disrupt ~ wCo/ft 


OC 


7WLr ’ 


(D14) 


which is the same as Eq. droli . Interestingly, Eqs. (ESI) 
and (1D14|) suggest that lower quench rates could be used 
when using high-viscosity fluids (larger rf). It has been 
reported that the binary liquid nitromethane-ethanediol 
is more forgiving in bijel fabrication than the W-L system 
(the viscosity of ethanediol is 16 times larger than for 
water) [26j . 


where 77 is viscosity, A r the typical length scale of 
the disruption (independent of particle radius) and v ~ 
A/At d . Combining Eqs. (|D 6 [) and (1D7[L we get 

At d oc -L - , (D8) 

7 WL r 

which is Eq. (|5]) . 

Alternatively, consider the equation of motion of a 
damped oscillator (compare Eq. (jD3[) h 

mil + fjbu-\- 4 ku = 0 , (D9) 

in which /i is a drag coefficient. We assume here that, 
at least initially, the drag mainly comes from the bulk 
fluids. In that case, 

H = 77 A 3 . (DIO) 


2. Jamming Time 

Consider a bijel surface S of area A(t), i.e. the area 
of the liquid-liquid interface between the two channels 
is decreasing during coarsening. Then the 2D packing 
fraction of particles on S is 

m = ’ (D15) 

with < 2 wl($) the particle-interface cross-sectional area 
and N the number of interfacial particles. Here, we as¬ 
sume that both clwl{0) and N are constant during the 
crucial (jamming) stages of bijel formation, for there is 
hardly any area left on S for new particles to attach to. 
Eq. (jD15[) still holds for the bijel in its final i.e. jammed 
state, so 


In our experiments Lf >> r, but if Lf ^ r then bulk drag 
may no longer dominate and effects of surface viscosity 
would have to be considered (which is outside of the scope 
of the current paper). 


4> f 


VawiA#) 

At 


Na WL (0) 0fA f , 


(Die) 
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which leads to 


= <t>t 


At 

A(t ) ‘ 


(D17) 


As it is L(t) rather than A(t) that is typically reported 
from simulations and experiments, we write 


A{t) — c g 


V c 

m ’ 


(D18) 


If 0i n is the packing fraction at which interfacial par¬ 
ticles start interacting, thereby affecting the phase sepa¬ 
ration [38|, then 

0f - 0in = (xi) /t in 

~ (lt) 1 11 * * * Vl ( t{ “ iin ) (D20) 

Atj = t{ - t in « (l - ^f) (lr) > 


in which c g is a geometrical pre-factor and 14 is the total 
volume of the bijel channel (which is constant during the 
phase separation of a symmetric binary liquid). Combin¬ 
ing Eq. (jD16[) with (jD18[) gives 


4>(t) 


i c g V c L(t ) 
L f CgV c 



(D19) 


where we have assumed that c g is constant i.e. the topol¬ 
ogy of the bijel does not change substantially during the 
final stages of (successful) formation. 


where in the second line we have used L(t) oc £, which 
is valid in the relevant phase-separation regime for bijel 
formation (viscous-hydrodynamic) (28) . 

Note that Eq. (jP20[) can explain several observations. 
First, the larger Lf, the longer the jamming time, which 
may help explain the empirical upper limit to bijel chan¬ 
nel width [23 • Secondly, the larger the coarsening speed 
vl , the shorter the jamming time. As vl increases with 
heating rate, through its dependence on the temperature- 
dependent interfacial tension [39], this may help explain 
why heating faster facilitates successful bijel formation 
(even for MPs). 
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